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Abstract 

A Markovian model for a quantum automata, i.e. an open quantum dynamical 
system with input and output channels and a feedback is described. A multi-stage 
version of the theory of quantum measurement and statistical decisions applied to 
the optimal control problem for quantum dynamical discrete-time objects is de- 
veloped. Quantum analogies of Stratonovich non-stationary filtering and Bellman 
quantum dynamical programming for the time being discrete are obtained. 

The Gaussian case of quantum one-dimensional linear Markovian dynamical 
system with a quantum linear transmission line is studied. The optimal quantum 
multi-stage decision rule consisting of the classical linear optimal control strategy 
and quantum optimal filtering procedure is found. The latter contains the opti- 
mal quantum coherent measurement on the output of the line and the recursive 
processing by Kalman-Busy filter. 

All the results are illustrated by an example of the optimal control problem for 
a quantum open oscillator at the input of a quantum wave transmission line. 



1 Introduction 

High perspective of applying quantum coherent electromagnetic generators of optical and 
infra-red frequency band for communication and control of quantum dynamical objects 
stimulates an increase of the interest in theoretical investigations of potential possibilities 
of information systems containing quantum channels. 

Due to fundamental limitations of quantum-mechanical measurement a specific prob- 
lem of optimal nondemolition measurement on the input and the output of quantum 
channels arises in such investigations. Here we shall consider such a problem for the 
channels with a feedback, corresponding to the optimal control in quantum open systems. 
It is essential in quantum theory that systems under the observation should be open, i.e. 
matched with channels, in order not to demolish them, by letting out an information. 



This paper gives the positive answers in a mathematically constructive way to the fol- 
lowing fundamental questions of quantum systems theory: Is it possible at all to observe 
and control a quantum dynamical system in the real time without not destroying it? If 
yes, what are the optimal strategies of that observation and control? How the dynamics 
of a quantum system is to be changed under the obtained information and its use as a 
feedback? What is are the fundamental limitations of quantum observability and control- 
lability? Is there any possibility to obtain a time continuous limit of such observation and 
control in a quantum system? 

The non-dynamical problem of quantum measurement optimization formulated pri- 
marily for detection and estimation in the static quantum communication systems by 
K. Helstrom [1] was studied intensively by several authors [2-7] within the framework 
of the single-stage (static) quantum-statistical decision theory. The dynamic problem of 
quantum nondemolition measurement for communication and control has been studied in 
details by author even in continuous time [8,9] since the pioneering paper [10]. However 
the earlier paper [11] on the solution of the discrete-time problem of optimal measure- 
ment has never been published in full, in spite of the practical importance of this case 
for the digital communication and control in quantum channels with a feedback. The 
novelty of this paper was such, that only a few people working in the newly open area 
of quantum stochastic processes could appreciate it at that time, and it was too far yet 
from applications. Recently, however, in view of the new possibilities of quantum compu- 
tations, the interest to quantum theory of communication and control has been renewed, 
and the time development of the discrete models of quantum open systems for commu- 
nication and control became actual. Moreover, after the development of time-continuous 
theory of quantum nondemolition measurement and filtering within the quantum stochas- 
tic calculus approach [12], these models can be considered as discrete time analogous and 
approximations of this theory. The discrete time case is mathematically simpler, it doesn't 
need the theory of quantum stochastic integration, and might be considered on its own as 
a dynamical programming for quantum computations, or multi-stage variant of optimal 
quantum-statistical decision theory. 

The quantum dynamical programming for multi-stage optimal measurement problem 
can be considerably simplified due to assumption that not only the processing of the 
measurement results but also the quantum measurement itself may depend on all previous 
measurement results. It corresponds to the assumption that we can choose a quantum 
measurement apparatus on the basis of the previous measurement data separately at 
every instant in time. Though in reality it is possible to imagine such a situation only 
for a finite number of stages and a finite set of measurement results (time and measuring 
scale being discrete) , this extension of admissible measurement and decision procedures is 
mathematically very convenient and from the physical point of view is not contradictory. 
The choice of the measuring apparatus and of the observed data processing according to 
all previous measurement results on the whole defines the strategy in multi-stage quantum 
decision theory described here. Within the framework of such an approach the problem 
of quantum filtering of random signal sequences was reduced in [13] to the well-studied 
problem of the static optimal quantum measurement on every fixed stage with conditional 
a priori distribution depending on the previous observed data. 

Here we describe the multi-stage quantum statistical decision theory applied to the 



problem of optimal control of a quantum Markovian discrete time system with a matched 
quantum channel. This theory may be considered as a quantum (operational) analogue 
of the stochastic control theory, based on Stratonovich theory of conditional Markovian 
processes [14], and Bellman dynamic programming [15]. 

The optimal filtering and the control strategy are found here in case of one-dimensional 
quantum linear Markovian system with quantum Gaussian noises and the mean-square 
loss function both in the discrete and continuous time. 

In order to pose the problem of measurement and control correctly from the physical 
point of view , let us consider the following motivating example. 



2 Controlled quantum open oscillator with quantum 
transmission line 

We are going to give a Markovian model of the simplest quantum system with a com- 
munication channel: the quantum open oscillator matched with a transmission line. It is 
an excellent mathematical model of a single-mode antenna for quantum radiophysics and 
optical control and communication. 

Let x be an operator of complex amplitude of a quantum oscillator with Hamilto- 
nian flx*x, which satisfies the canonical commutation relations with x* being an adjoint 
operator 

[x, x*} = xx* — x*x = hi (2.1) 

where 1 is the unit operator, and h > is the Planck constant. 

Assume that in general case this oscillator is controlled by the complex amplitude u 
by means of a quantum- mechanical transmission line with wave resistance 7/2, where the 
operator of the wave y (t — travelling from the oscillator into the line is measured. 
In the simplest case of ideal conjugation between the line and the measuring apparatus, 
when there is no reflection of the wave travelling from the oscillator, i.e. in case of the 
matched line, x (t) and y (t) are described by the pair of linear equations [16] 

dx (t) /dt + ax(t) = 7«(i) + v(t), x (0) = x, (2.2) 

y (t) = ax (t) - dx (t) /dt = 7 (a; (t) - u (t)) - v (t) , (2.3) 

where, generally speaking, a is a complex number with fixed real part, a + a = 7, and 
with arbitrary imaginary part depending on the choice of the representation, v (t + 
the amplitude operator of the wave travelling out of line towards the oscillator, this 
operator is responsible for the commutator preservation. Under natural for super-high and 
optical frequencies assumption of narrowness of the frequency band which we deal with 
the commutators for v (t) in the representation of "rotating waves" have delta-function 
form [17]: 

[ v (t) , v (t')] = 0, \v (t) , v (If)*] = ^H15 (t - If) . (2.4) 



Integrating equation (2.2) and taking into account that v (t) does not depend on 
x(t') when t > t', it is easy to verify that the commutator [x(t),x(t)*] is constant, 
moreover, x (t) commutes both with y(t') and y{t')* when t > t', and the commutators 
for y (t) ,y (£') , y (t 1 )* coincide with (2.4). The latter means that considering von Neumann 
reduction which appears as a result of some quantum measurement of y (t) at previous 
instants of time t' < t does not affect the future behaviour of x (t l ) , y (t 1 ) , t 1 > t, so 
that equations (2.2), (2.3) remain unchanged. This fact together with the Markovianity 
hypothesis of the quantum process x (t) which hold for quantum thermal equilibrium 
states of the wave v(t) in case of narrow band approximation [17] simplifies to the large 
extent optimal measurement and control problems for the simplest quantum dynamical 
system mentioned above. 

Let us assume, that the initial state x is Gaussian with the mathematical expectation 
(x) = z and 

((x — z) (x — z)) = 0, ((x — z)* (x — z)) = hT>, 
v (t) is the quantum Gaussian white noise, which is described by the following correlations 

( v (t) v {?)) = 0, (v (t)* v (*')) = haS (t - t') , 

with <T = 7 (exp (TiVt/kT) — l)~ 1 for the equilibrium state with the temperature T, where 
k > is the Boltzmann constant. 

As an example, let us try to choose the optimal measurement of the controlled quantum 
oscillator (2.2) with transmission line (2.3), so that to minimize its energy Q (x* (t) x (t)) 
at the final instant of time t = r by means of the control strategy the norm J T \u(t) \ 2 dt 
of which should not be too great. If we want also to force the quantum amplitude x(t) to 
follow the classical process u(t), this problem can be characterized by the quality criterion 

ft (x (t)* x (r)) + I" (6u (t)* u(t)+u (x (t) - u (t))* (x (t) - u (*))) dt. (2.5) 
Jo 

Here 9, uj > are parameters responsible for the measurement quality: when 6 = ft = 
(2.5) corresponds to the problem of pure filtration, when u = 0, 9 ^ 0, it corresponds to 
the pure control problem. 

It will be shown below (see §5) that the optimal measurement minimising criterion 
(2.5) is statistically equivalent to the measurement of the stochastic process z (t) — x (i) + 
x° (t) described by Kalman-Bucy filter: 

dx (t) /dt + ax (t) = 7M (t) + k (t) (y (t) - 7 (x (t) - u (t))) . (2.6) 

Here x (0) = z, k (t) = (7S (t) — a) / (/j, + a) , S (t) is the solution of the equation 

d£ (t) /dt = (a - 7 E (*)) (ji + 7 S (*)) /(ji + <t), E (0) = S, 

dx°(t)/dt + ax°(t) = K(t)(v°(t) - 7 x°(t)), x° (0) = 0, (2.7) 
where v° (t) is the amplitude operator with commutators 

[v° (t) , v° (t')\ = 0, \v° (t) , v° (if)*] = -h^5(t - t') 



which change the quantum process x (t) into the classical (commutative) diffusion complex 
process, and with correlations of vacuum noise of the intensity ji = if 7 < and /i = 7 
if 7 > 0: 

(v° (t) v° (0) = 0, (v° (t)* v° (0) = hfi5 (t - f) . (2.8) 

For instance, such measurement takes place by the heterodyning [7] where v (t) stands 
for a standard wave. In this case the optimal control strategy u° (t) coincides with the 
classical one: u° (t) = — A (t) z(t), where \(t) = ( 7 Q (t) — u) / (9 + u), and Vl (t) is a 
solution of the equation: 

-dtt{t)/dt = (w- 7 fi(t)) (e + jn(t))/(e + u), n(r) = n, (2.9) 

which together with (2.7) defines the minimum quantity of losses (2.5): 

n (n (o) e + £ (n (t) a + ( 7 n (t) - uf s (t) /{e + u)) dtj + n (o) | z \ 2 . 

By setting a = 0, u; = 0, we obtain in particular the solution of the terminal control 
problem for an oscillator with thermal noise equal to zero. But in this case unlike the 
classical one the optimal measurement remains indirect and the equation (2.7) remains 
regular corresponding to the white noise in the channel of intensity | 7 | h. Thus to 
consider the quantum measurement postulates is statistically equivalent to the adding of 
white noise into the channel of intensity | 7 | h what excludes the singular case of pure 
measurement of the amplitude x. 

It is interesting to note that in the case of thermal equilibrium when 7 > 0, T > and 
£ = (exp {TiVt/kT} — 1) _1 the optimal amplification coefficient k (t) equals to zero which 
means the possibility of optimal control of the quantum oscillator without measurement. 
It also holds when u = 7 fi, the solution of equation (2.9) is stationary and optimal feed- 
back coefficient A (t) equals to zero. But in the contrary case 7 < 0, T < which corre- 
sponds to the active medium of the oscillator (laser) the optimal coefficients k (t) , A (t) are 
strictly negative and non-zero even for the stationary solution S (t) — 0, Q (t) — 9/ | 7 | 
of equations (2.7), (2.9). 

3 Quantum dynamical filtering 

Now let us give a rigorous setting of the quantum dynamical observation problem for the 
optimal control of a quantum-mechanical object when time is discrete t G {t k } k=o,i,...- 
Let Ak be von Neumann algebras on a Hilbert space H, each is generated by one or a 
few dynamical variables (operators) Xk = {.x l k ) %e in H. One can consider a quantum- 
mechanical object in the Heisenberg picture at the instant of time tk+i > t k > with 
Xk = x(tk), such that all algebras Ak are equivalent to the initial algebra Ao = A, 
generated by the positions and momentums x = (q,p) at t = 0. Let Bk,k = 1,2,... 
be von Neumann algebras of observables generated in H by output dynamical variables 

Uk = {ykf ' by means of which this object can be observed in a nondemolition way 
say, on the time intervals (tk-i, tk]- As it has been shown above on the example of the 
matched transmission line, the output observables bk G Bk in the matched channels should 



commute with all present and future operators a k i G , k 1 > k of the dynamical system, 
but not necessarily with the past ones a k > G A k >, k < k. This commutativity condition 
together with the commutativity b' k b k i = b k ib' k for all b' k G B k , b k < G B k > \fk' ^ k will be 
referred as the nondemolition condition. 

Let us denote V k , TZ k the dual spaces to A k , B k with respect to some standard pairings 
< ., . >, say the subspaces of trace class operators n k G A k , p k G B k which are dual to the 
simple algebras of all bounded operators a k G A k , b k G B k on the corresponding Hilbert 
spaces with respect to the bilinear trace-forms 

< 7i k , a k >= tr [n k a k ] , < p k , b k >= tr [p k b k ] , 

and denote S k the corresponding subspace dual to the von Neumann algebra B k V A k 
generated by the commutating B k and A k . We shall use the operational terminology, 
briefly summarised in the Appendix. Thus we shall call the positive normalized elements 
K k G V k , p k G TZ k and a k G S k , which are usually described by the statistical density 
operators, the statistical states of the quantum object at the instants of time t k , the states 
of the channel on the interval (t k -i, t k ], and the joint state of the object and channel at 
the moment t k respectively, or simply the states on A k , B k and B k V A k C B k ® A k . 

Now we adopt the hypothesis of Markovianity of the Heisenberg dynamics, restricted 
to the described quantum object and output channel in TC, with respect to a given state 
of the whole system uj. Let all the induced states a k = cu\ {B k V A k ) , k = 1, 2 . . . and 
their restrictions p k , n k on B k , A k be defined by the initial state ttq = tt on A® = A and 
by a family {M k } ^-1,2... of statistical morphisms ir k ^i 1— > o k = -n k -\M k . These transition 
maps V k -\ — > S k can be described as the pre-dual to positive normalized superoperators 
M k : B k ® A k — > Ak-i having for the simple algebras the form 

M k c k = tr B o [p° k Ck ] , Vc fe eB k ®A k . 

Here p k , k = 1, 2, . . . are states on some algebras B k , for which the simple algebras B k ®A k 
are isomorphic to the von Neumann tensor products A k -i ® B° k , and tr^o is the partial 
trace on B° such that M k [a fe _i ® b° k \ =< p° k , b° k > a k ^ for all a fc _i G A k -i, b° k G B° k . This 
assumption corresponds to the requirement that the channel should be matched with 
the object and implies the semigroup dynamics [20] n k -i 1— > n k = trg k {n k M k } of the 
quantum-mechanical object with discrete time. Furthermore, we shall suppose that every 
morphism M k may depend on the results ( k = {(k'} k > <k of previous measurement data 
( k r G Z,k' < k, say via dependence of some controlled parameters u G U of the sequence 
{Ck'} k / <k due to a feedback ( k 1— > u. 

The nondemolition measurements during the time intervals (tk-i, t k ] are described by 
positive operator-valued measures b k (d() G B k , k = 1,2,... on the data space Z 3 ( 
with a given Borel structure of the measurable subsets dz C Z such that b k (Z) = 1 
is the identity operator of B k . We shall assume that every Z-measurement b k (d() also 
may depend on all preceding measurement results (1, . . . , ( k -i, and not only due to a 
dependance onti G f/ and the feedback, but directly, being adaptive in time. The functions 
( k 1— > (Mfc (( k ^J , b k (( k ,d(j) are supposed to be weakly measurable in the sense that for 
all 7r fe _! G Vk-i and a k G A k and all Borel subsets d( C Z the complex functions 



are Borel functions on Z k = Yl k / <k Z k r, where Z k = Z,Zq = U. We shall call every 
sequence {b k (C fc )^C)} fc 12 °f suc h "conditional", or adaptive measurements the mea- 
surement strategy. 

Let us denote B k (( k ,d(j the conditional transition measures Vk-i — ► V k , that is the 
operational-valued conditional measures on Z, defined as the predual to superoperator 
values Bk (( h , d(j : A k — ► Ak-i by the formula 

a k » B k (C fc , d() a k = M k (C fc ) [h (( k , d() ® a fc ] , (3.1) 

and denote ir k (d( k+1 ^j the V k - valued measures on Z k+1 obtained for k — 1, 2, . . . by the 
recurrency 

7T fc (d( k x dC) = TTfc— i 5fe (C", dC) (3-2) 
with the initial condition 7r (^C 1 ) — n $ ( u o, dC 1 ) if Z 1 = U. 

Lemma 1 All the measures n k (d( k+1 ^j are positive in the sense that 

J<7c k (d( k+1 ) , a k (C fe+1 ) > > 

for all A k - valued positive measurable functions a k > and are normalized, 

< n k (Z k+1 ) , 1 >= 1, where 1 is the identity operator of Ak- 

Proof As the superoperator- valued measures B k (( k , d(j are positive and normalized in 

the sense that J B k (( k , d() [a k (( k ,()] > for all a k (C fc+1 ) > and B k (( k ,z) 1 = 1, 
the lemma can be easily proved by induction, using the positivity and normalization of 
7r . Thus the measure 7r k (d( k+1 ^j, obtained by the recurrency (3.2), describes the total 

statistical state on the algebra A k and on the expanding space Z k+1 = Z k x Z ■ 

Let us define a posteriori state of the object at time t k as V k -valued Radon-Nikodim 
derivative 

TTfc-l (f) = 7r fe _! (d( k ) I < T^ k —\ (d( k ) , 1 > (3.3) 

which exists in the weak sense due to absolute continuity of it k -i with respect to < 

7Tfc-l, 1 > 

Theorem 2 The a posteriori states n k (C fc+1 ) , k — 1,2, . . . can be obtained by the non- 
linear recurrency 

*k (C*, C) = 7r*_i (C fc ) T k (( k , C, Tik-i (C*)) , (C 1 ) = 7T, 
where T k ((, k+1 , ~K k -i) is the (V k -i — > V k )-valued Radon-Nikodim derivative 
T k (C k , C, ^-i) = B k (C h , d() I < Kk-iBk (C h , d(),l>. 



Proof The nonlinear transition operations T k are defined in the weak sense almost ev- 
erywhere by the Radon-Nikodim derivatives 

< Kk-iTk (C k+ \ n k -i) , a k >=< 7r fc _i£ fc (( k , d(),a k > / < n k -iB k (( k , d() , 1 > . 

The proof of the theorem follows immediately by induction due to the Bayes formula 

< 7T k (d( k X d() , 1 > / < 7Tfc_i , 1 > = < TTfc-! (C fc ) 5fe (C", dC) , 1 >, 

from the definitions (3.2), (3.3) ■ 

Note that the equation (3.4), describing the conditional Markovian evolution of a pos- 
teriori state of a quantum-mechanical object, can be regarded as a quantum generalization 
of Stratonovich nonlinear filter equation with discrete time. A semi-quantum case when 
a partially observed object is described by a classical Markovian process {x k } k=0 1 and 
the channel is non-classical, was considered in [12]. 



4 Quantum dynamical programming 

Let us consider the problem of optimization of the observation strategy \b k (( k , ^C)} on 
the fixed discrete time interval [0,-fT]. The optimal strategy {K} k e[oK) * s defined as a 
strategy, which minimizes the average cost 



K 

a =< tt k , a K > + J2 / < Tfc-i (d( k ) , Cfc-i (C*) >, (4.1) 
k=i J 



given by a self- adjoint semi-bounded operator clk £ Ak of final losses, and by similar 
operator-valued functions ( h+1 h- > c fe (C fe+1 ) £ ^4fc, k — 0, . . . , K — 1. (In the case of 



unbounded and on ly their spectral measures should belong to Ak and A k .) 

Let us remark that the cost (4.1) does not depend on the last measurement b K (C, K , d(j 
which can be chosen arbitrarily, and t\k — [Z k+1 ). As it follows from definitions (3.1), 



(3.2) the YX'=i m (4-1) for any k — 1, . . . , K is independent of the measures b k i 

for k 1 > k. Hence in order to find the optimal Z-measurement b k (( k , d(j from some 
k < K it is enough to vary the future average observation cost functional 

K f 

a k =< ttk, a K > + J2 < ^k'-i (d( k ') , cv_i (C*) > • (4.2) 

k'=k+l 

Lemma 3 The explicit dependence of a k on b k (( k , d(j is affine 

«*=/./ < Pk (d( k , C) , b k (( k , d() >, (4.3) 



where p k (d( k , () = vr fc _i (d<C fc ) A k {( k , () ■ Here A k (( k+1 ) is a (V k -i -> H k ) -valued func- 
tion on ( k+1 which is defined as predual to the superoperators 

b k ^A k (( k+1 )b k = M k (( k ){b k ®a k (( k+1 )}, Vb k eB k , (4.4) 



where a k (C fc+1 ) is an operator-valued function on Z k satisfying the linear inverse-time 
recurrency 

o fc -i (C*) = / B k (C,\ dC) a k (C fc , C) + Pik-i (C") , (4.5) 
/c = 1, ...,K with the boundary condition ax (C^" 1 " 1 ) = clr- 
Proof First let us prove that the future losses (4.2) can be represented as 

ttfc = J J < TTfc (d( k x d() , a k (( k , () >, 

where a k (C fc+1 ) £ A k is the solution to the equation (4.5). It is obviously valid for k — K, 
and if it is true for a k < K, then substituting (3.2) into this representation of a k , we 
obtain 

J <n k (d( k+1 ) , a k (C fc+1 ) > + J < TTfc— i (d( k ) , c k - X (( k ) > 
= J zk < 7r fc _i (d( fe ) , £ 5fe (C*, dC) a* (C, C) > +c*_i (C") • 

So this is also valid for a. k _\ with a^-i given in (4.5), and by using the inverse-time 
induction, it is valid for any k G [0, K). Now we can obtain (4.3) by 

< n k (d( k x d() , a k (C fc , C) >=< «k-i (d( k ) B k (c\ d() , a k (( k , () > 
=< 7r fc _! (d( k ) M k (( k ) , b k (C fc , dC) ® afc (C k , C) > 

=< vr^ (d( k ) A k (C k ) , 6 fe (C fe , d() >=< Pk (d( k , C) , b k (C k , d() >, 

where we used the definitions (3.1) and (4.4) for the operations B k and A k ■ 

Theorem 4 If the strategy \b° k {(, k , d(j } feg j 1 ^ ^ s optimal for the cost functional (4-1), it 
satisfies the following system of equations 

(p k (d( k , C) - A fc (d( k )) b° k (C fe , dC) = 0, (4.6) 

fc G [1, if), where 

h(d( k ) = I Pk(dC k ,()b° k (C k ,d(). 

These equations together with the system of inequalities 

Pk(d( k ,()>\ k (d( k ), ke[l,K) (4.7) 

give the necessary and sufficient conditions of the optimality for quantum measurement 
strategy b° k , k = 1, ...K — 1 corresponding to the minimal values 

a° k = J <X k (d( k ) , 1 > (4.8) 

of the future average costs (4-2). 



Proof As the variables b k \( k , d(j , k — 1,2, . . . , K oi the functional (4.2) are indepen- 
dent, the optimal measure b° k {(^ k , d(j minimizes the affine functional separately for every 
fixed family (C fc \^C)} fcl fc - The necessary and sufficient conditions (4.6), (4.7) of 

optimality for b° k (( k , d(^j, minimizing the affine functional (4.3) with a fixed k, follow im- 
mediately by the linear programming method, as it was noted in the single-stage theory 
of optimal quantum measurements [2, 4-7] ■ 

Note that the minimal value a° of the total average cost (4.1) is given by the solution 
a° = of the recurrency (4.5) with B k = B k at k = as a° =< n, a° >. 

Let us note that with the help of the a posteriori states n k (C fe ) > one can write condi- 
tions (4.6), (4.7) in the following form 

{Pk(( k ,()-\k(( k ))b° k (( k ,d()=0, (4.9) 

Pk (C k+1 ) > X k (C fc ) , ke[l,K), (4.10) 

where p k ((, k+1 ^j = Kk-i (C fc ) A k (( k+1 ^J. According to the Bellman dynamical program- 
ming method [15] the verification of the optimality condition formulated above can be 
carried out sequentially in inverse time k = K — 1, . . . , 1 applying the recurrence (4.5) for 
the superoperator A k ((, k+1 ^j after solving the filtering recurrent equation (3.4). 

The optimal control of Markovian partially observed quantum-mechanical object can 
be reduced to the optimal measurement problem investigated above as follows. Let 
M k (wfc_i) : Vk-i — > T^k <8> Vk be the quantum statistical morphisms (transitions) con- 
trolled by some parameters u k G U, k = 0, . . . , K — 1. A control strategy {lk} k< K is given 
by a choice of the feedback, defined by a measurable dependence "f k of each u k on all mea- 
surement data r\k> k' < k, and also on the preceding controls uy, k' < k. The optimal 
control for a fixed measurement strategy is supposed to minimize the average cost defined 
by a final operator ax and operator- valued cost functions Ck (uk) , k — 0, . . . , K — 1. De- 
noting C 1 = uq, ( k = (uo,r)i, Mi, ... , rjk-i, Uk-i), ( = (v, u ) ? the average cost functional 
even with random control strategies can be represented in the form (4.1), given by the 
quantum measurement strategy \bk (( k , ^C)} on Z = Y x U of the form 

b° k (( k , drj x du) = b° k (C fc , drj) 5 ( 7fe ° (C fc , v) , du) (4.11) 

and c (C 1 ) = c (u ), c k ^C fe+1 ) = c k (u k ). The quantum optimal control problem can be 
formulated then as one of searching for the optimal Y x U -measurements b° k {(, k , d(j , k G 
(1, k) , and an optimal initial control u° corresponding to the minimal value 

a° — inf < Ai (u) , 1 > + < 7r , c (u) > 

of average cost (4.1). In general, the optimal measurement strategy may not be in the 
product form (4.11), but if there exists a non-randomized strategy u° k = ^° k (( k , i]j , k G 

[1, K) for some F-measurements b° k (C, k , drfj for which the Y x U - measurements are 
optimal, where S(-,-) is the Dirac 5- measure, then the data spaces Y may be called 



the sufficient spaces. The optimal measurements b° k [( k , drjj on sufficient data spaces Y 
satisfy obviously the equations 

[pk (C", V, ll (C*, r,)) ~ A* (C*)) b° k (( k , dri) = 0, k G [1, X), 

where 

A* (C*) = /p* (C\ V, ll (C*, vj) K (C k , dri) , 

which together with the inequalities (4.10) are necessary and sufficient for the non- 
randomized control strategy {~i k }- 

5 Quantum filtering in Boson linear Markovian sys- 
tem in a Gaussian state 

We examine a Markovian one-dimensional quantum dynamical system, described at dis- 
crete instants t k = kA by the algebras At and B k , which are generated by the non- 
selfadjoint operators x k 7^ x k and y k 7^ y k respectively, satisfying the canonical commu- 
tation relations. Let us suppose that they act in the same Hilbert space 7i, where they 
satisfy the linear quantum stochastic equations 

x k = + f3u k _i + v k (5.1) 

y k = 7x fc _i + 5u k ^ + w k . (5.2) 

Here 0, f3, 7, 5 are some, in general complex parameters, the controls u k can also accept 
complex values, xq = x is the initial operator in 7i, generating the algebra A, and Vk,w k 
are some operators in 7i, generating the algebras B k . To obtain the Markov dynamics, 
we need to assume the independence of x and all the pairs (v k , w k ) such that the algebras 
A and B k , corresponding to different instants of time tk, commutate, and the joint state 
uj is the product of the states on A and all B k , k = 1,2.... We shall define the canon- 
ical commutation relations for the generating operators x,Vk,Wk with their adjoints as 
following: 

[x,x*] = hi [v k ,v* k ] = (i-\ ( f ) \ 2 )hi, 

[ w k,wl] = (e- I 7 I 2 ) KX] = -07^1, (5-3) 

where h > and 1 is the identity in H ( other, unwritten commutators, including all 
those corresponding to different instants of time to be equal to zero.) Here the choice of 
the commutator [wk,vl], responsible for the commutativity [yk,x* k ] = is essential, the 
other nonzero commutators are chosen so that the commutators 

[x k ,x* k ]=hl, [y k ,y* k }=ehl 

should be constant. The described system we shall call the discrete linear Markovian 
quantum open oscillator. 

Let us describe the states ir k G Vk by the Glauber [21] distributions pk (£), £ G C, nor- 
malized on the complex plane C with respect to the Lebesgue measure o?£ = o?Re£ cflm£ / irh. 



In the representation described in the Appendix, the Markovian morphisms V k -i — > V k , 
corresponding to the linear equations (5.1), (5.2), transform the distributions pk-i (£) into 
the two-dimensional distributions 

9k (£, V)=jqk(t- ^ -Pu,ri- ^ - 5u) (e) d£\ (5.4) 

where q k (£,rj) are some other (not necessarily Glauber) distributions on C 2 , which de- 
scribe the independent states p° k on algebras B° k . 

When e — 0, the operators y k ,y k are simultaneously measurable, and the a posteriori 
states on A k under the fixed spectral values y k > = r\y and u k ', k' < k are defined recurrently 
by the a posteriori Glauber distribution p k | C fc_1 ) according the Bayes formula 



Vk 



Here g k (^,rj | C fc )are the distributions obtained by substitution of \ ( k ^j into (5.4) 

instead of p k -i (£), and 

r k (r 1 \C k )= j 9k (t,v\C k ) 

are the probability distributions, describing the complex values r] k , which arise as the 
results of the direct measurements of y k under the fixed ( k = (u , rji, ui, . . . , r) k -i, u k -i). 

When £ ^ 0, only indirect measurement of y k are possible which are described, for 
instance, by the B k -valued measures 

h (drj) = #m fe (77 - y k ) #drj, (5.5) 

represented by some distributions m k {rj) on C as it is described in the Appendix (A. 3). In 
this case in order to calculate a posteriori Glauber distribution one should change q k (£, rj) 
in formula (5.4) for the distribution 

Qk (f , V) = J m k (v - r] 1 ) q k (f , V 1 ) dr] 1 . (5.6) 

Theorem 5 Let the initial state n of the quantum oscillator be described by the Glauber 
distribution of Gaussian type 

p(0 = ^exp{-|e-^| 2 Ms}, (5.7) 
the distributions q k describing the transitions (5.4), be also Gaussian: 



q k (£, r,) = v-.s. v, .-m,/^ (58) 

ov— I V \ z 

and the measures b k are described as in (5.5), by the Gaussian distributions 

(77) = — exp |— I 77 | 2 /h/i} . (5.9) 



Then a posteriori states (3.3) at each instant k = 1,2, ... , are given by the conditional 
Glauber distributions of Gaussian form 



Pk U I C fe+1 ) = =- exp {- | £ - | 2 /hY, k \ , (5.10) 

^k 

where z k , T, k are defined by the recurrent equations of the complex Kalman filter: 

z k = 0z fe _i + + K k (r} k - izk-x - <Jttfc_i) , z = z, (5.11) 

£ fc =| I 2 s fc-i + I «* I 2 S = S, (5.12) 

where 

K k = (07S fc _i - v) ^ k =| 7 | 2 + I/ 1 , V 1 = V + fJ > . 

Proof Due to the chosen representation, the proof is similar to the derivation of the 
classical one-dimensional Kalman filter for the complex Gaussian process x k given by 
(5.1) and y k = y k + w%, where wl are independent Gaussian variables with zero mean 
values and the covariances /j, > e. (For this proof see, for instance, [22].) One should only 
take into account that distributions (5.6) are also Gaussian of the type (5.8) with the 
parameter z/ 1 = v + \x instead of v. Substituting q(£,,i]) in (5.4) by q 1 [C,,i]) and p k -\ (£) 
by the conditional distribution p k -\ \ C fe ) of type (5.10), we obtain 

where p k (j; \ ( k l ^j is the Gaussian distribution (5.10) with the parameters (5.11), (5.12), 
and 

rl (p | ( k ) = ^- exp {- | 77 - 7^-1 | 2 /h* k } . (5.13) 

^ k 

Thus the quantum Gaussian filtering is controlled by the classical Kalman filter for the 
complex amplitude in the Glauber representation ■ 

Note, that in distinction from the classical case, the covariance matrix of distributions 
(5.8), (5.9) should not only be non-negative definite but should also satisfy the Heisenberg 
uncertainty principle 

as it follows from inequality (A. 5). In particular it excludes the case \i = of the direct 
observation of y k when e > 0. 

As shown in the next paragraph, a posteriori mathematical expectations z k with fi = 
max (0, e) appear to be the optimal estimates u k = z k of the operators x k with respect 
to the square quality criterion c k (u k ) =: \ x k — u k | 2 : with the minimal mean square error 
foSfc. In the commutative l\ — this optimality was proved in [11]. 

Note, that instead of calculating z k by means of the recurrent formula (5.11) using the 
results (rji, ...,r) k ) of the indirect measurement (5.5) one may regard z k itself as a results 
of the measurement described by the B k -valued measure: 



h (( k , dzj = #n k (z - x k ) #dz, 



(5.15) 



where 

Uk ( Z ) = I \2 mk ( Z / I Kk I) ' 

I K k | 

and 

x k = (pZk-x + (3u k ^x + n k (y k - -fz k -i) (5.16) 

is an operator, depending on the values z k -i,u k , and independent of the preceding mea- 
surement and control results. 

It is interesting to consider the time continuous limit, when the quantum oscillator 
(5.1), (5.2) is described by the quantum stochastic differential equations 

dx (t) + ax (t) dt = j3u (t) dt + v (dt) , (5.17) 

y (dt) = 7a; (f ) dt + 5u (t) dt + w (dt) , (5. 18) 

i.e. by equations (5.1), (5.2) with x(t k ) = x k , y (At k ) = y k , 4> — 1 — «A, (5 ~ 
(3 A, 7 ~ 7A, e ~ eA, where (Atk) = t k — t k -i = A tends to zero. In addition to that 
the commutation relations (5.3) change in the following way 

[x, x*} = HI, [v (dt) , v {dt)*\ = (a + a) hdtl, 

[w (dt) , w (dt)*] = ehdtl, [w (dt) , v (dt)*] = -7/idtl, 

and the other commutators including those corresponding to the different instants of time 
are equal to zero. By passing to the limit as A — > when a ~ a A, v ~ v-A, v ~ u-A, 
it is easy to obtain under the assumptions of the Theorem 3 that a posteriori state it (t, £*) 
is described by the Glauber distribution p(t,£ | £*) of Gaussian type (5.10) with the 
parameters z (t) , E (t) which correspond to the Kalman-Busy filter 

dz (t) + az (t) dt — n (t) (rj (dt) - (72 (f) - 5u (t)) dt) . (5.19) 

Here re(t) = (7E (t) -v)/u l , u 1 = v + //, z(0) = z, E (0) = E, 

dE (t) /dt + (a + a) E (t) = a- \ x (t) | 2 v 1 , 

and 77 (rft) are the results of the corresponding indirect measurement of y (dt) which are 
realized by the measurement of the sum y (dt) + w° (dt), where w° (dt) is an independent 
quantum white noise, defined by the coefficients e, jj, : 

[w° (dt) , w° (dt)*] = -ehdtl, (w° (dt)* w° (dt)} = fihdt. 

As shown at the end of the next paragraph, such "continuous" measurement appears 
to be also optimal in the Gaussian case when \i = max (0, e). 

6 Optimal measurement and control for quantum open 
linear system 

In the following theorem it is not required that the distributions po, q k and m k should be 
Gaussian and it is assumed only that they should have the zero mathematical expectations, 
and the covariations should coinside with the covariances T,,S,v,u,ij, of the distributions 
(5.7) - (5.9) respectively, not necessary being of the form (5.11). 



Theorem 6 Let the operator of final losses be quadratic: ax = VLx* k Xk, where Vt > 0, 
and 

c k (u k ) = UJX* k X k - $u k x k - $u k x* k + ■d 1 | u k | 2 , u > 0, ft 1 > (6.1) 

be quadratic loss operators for all k e [0, K). Suppose u k = —\ k z k , k G [0, K) is a linear 
control strategy, where z k are the linear estimates (5.11) based on the results r\ k of the 
indirect measurement (5.5), and 

\ k = (Wt+i - ■&) /T fc , (6.2) 
with Tfc =| (5 | 2 fifc+i + ft 1 and Q k satisfying the following equation 

Q k =| | 2 fifc + i + to- | Afc | 2 Tfc, = Q. (6.3) 
TTien £/ie operators of future losses (4-5) are also quadratic: 

a-k (( k+1 ) = 41 + Tfc | Mfc + \ k z k | 2 +Q k x* k x k 
+ Tfc (z k - x k )* (z k - x k ) - 2ReAfc (u k + \ k z k )* (z k - x k ) , (6.4) 

where 

K 
i=k+l 

Tfc =| Afc | 2 Tfc+ | - /tfc +1 7 | 2 Tfc +1 , r fc = 0, (6.5) 

and 

Afc = <p(3tt k+1 - 

Proof In the representation 

a K =: a K (x K ) :, c k (u k ) =: o (x k , u k ) :, a k =: (a: fc , C fe+1 ) : 

the recurrent equation (4.6) has the form 

ttfc (6, C fe+1 ) = / afc+i (f , C fe+1 , »7, u) q 1 (f - 06 - V-lik- Su k ) d£dr] + a (f fc , u fc ) 

(6.6) 

where 

it = -Afc+iz, z = 0Zfc + /3«fc + /tfc+i (77 - 7^) . 

Let us assume that the function a k (£) has the quadratic form (6.4); in particular, it 
has this form at k = K, namely ax (£) = ^ | £ | 2 - Inserting the latter into (6.6) and 
integrating, we obtain, that the function a k -i is of the same form with Tfc_i = \ (5 | 2 
fifc and Qfc_i,r fc _! given by (6.3) and (6.5), and 

4-i = d k + h (fl k a + F k (a + 2ReK k v + v 1 \ n k | 2 )) . 

Summing Yif =k (4-i — 4) and taking into account that d K — and IV = 0, we obtain 
(6.4) also for k — 1 ■ 



Lemma 7 Let us assume that starting from the instant k + 1, the controls Uk are chosen 
to be linear u k i = —\ k iz k i with the coefficients (6.2), where z k i, k 1 > k depend linearly on 
the results of the subsequent indirect measurement rj k+ i, . . . , r/ K _i by virtue of the formula 
(5.11) with the initial condition z k = z. Let also the indirect measurements be described 
by the Gaussian distributions (5.9) up to the k. Then the operator p k (( k ,(^)> defined in 
(4-9), has the following normal form 

Pk (c fe ,c) = A fe (c k ) + 

+ : (r k \u + X k x k | 2 + I - Kfc+i7 | 2 T k+ i \z-x k | 2 ) r° k (y k | :, (6.7) 



where 



A fc (c k ) (n k \x k \ 2 + (h (n k + r fc ) s fc + d k ) 1) r° k (y k I c fe ) :, 



the operator x k , defined in (5.16), is linear with respect to y k , and r^{r] \ ( h ^j is the 
distribution (5.14) with the parameters v° = v + p°, where p° = max (0, e). 

Proof Indeed, the operator p k similar to the density operator p e lZ k is defined 

by the distribution r k (rj, C fc+1 ) = 

J j ' a k (f , C k+1 ) Qk (i - - I3u k -i, r] - 7^_i - 5u k -i) | C fc ) dtdZk-i. 

(6.8) 

It is a symbol of the contrary order (see the Appendix), which is normal when e < and 
antinormal when e > 0. In the former case, inserting the operator symbol (6.4) into (6.9) 
and integrating with respect to the Gaussian type of the distribution p k -i (t;k-\ I C fc ) > we 
obtain (6.7), where r° k (rj \ ( h ^j coincides with the distribution r k (rj \ C fe )of the Gaussian 
type (5.13) with the parameter v 1 = v. In the latter case e > 0, the normal symbol of the 
operator p k (^C fe+1 ) is obtained from (6.9) by means of the convolution of type (A. 2) with 
the distribution (5.9) with p = e, and in the result of the parameter v increases for e. In 
this case r° k (??|C fc ) is also the normal symbol of the conditional density operator p k (( k ^ 



on B 



k 



Theorem 8 Let the quantum oscillator (5.1), (5.2) be described by the Gaussian initial 
and transitional distributions of the Gaussian form (5.7), (5.8), and the quality criterion 
(4-2) be defined by the quadratic final and transitional operators ax = flx* K xx and c k (u k ) 
of form (6.1) respectively. Then the optimal strategy is linear: u k = —\ k z k , where \ k is 
defined by (6.2), and z k are optimal linear estimates (5.11) based on the results {i]i} i<k 
of the coherent measurements (5.5) which are described by the distributions (5.9) with the 
minimal value of the parameter p = p°. 

Proof We should verify the necessary and sufficient optimality conditions (4.10), (4.11) 
for the operator (6.7) and the mentioned above measurement at each instant k. As 
Yfc, Tfc > 0, and the density operator : r k (y k \ ( k ^j : is non-negative definite, the differences 
Pk (C fe+1 ) — A fe (C fc ) are non-negative definite operators as well. It remains to verify the 



equations (4.13) for the optimal strategy 7£ [( k , rjj = —XkZk of the coherent measurements 
(5.5) or, what is the same, of the measurements (5.15) with the Gaussian distributions 
n% (z), corresponding to the case \i = \i° . Inserting u = —XkZ into (6.7) and taking into 
account (6.5), we obtain 

Pk (C k ,V,l°k (C\v)) - X k (C*) =T k :\z-x k | 2 r k (y k \ ( k ) : . 

Thus, equations (4.13) with e > can be written in the form 

(z - x k ) #n° k (z - x k ) # = 0, 

and the adjoint ones can be written for e < also as 

#n°(z-x k )#(z-x k )=0. (6.10) 

The operators j^n° k (z — Xk)# described by the Gaussian distributions n° k (z), which re- 
alize the lower bound of the Heisenberg inequality, are well known as proportional to 
coherent projectors [8]. The operators x k when e > 0, are proportional to the annihila- 
tion operators, and when e < 0, they are proportional to the creation operators, for which 
the coherent projectors are the right and the left eigen-projectors respectively. Hence, 
the equations (6.9) is satisfied in the case e > 0, and the equation (6.10) is satisfied if 
e < 0. Note that, in the antinormal case when the coherent projectors are described by 
the Dirac distributions 5 (z) on C, these equations are written as the identities 

(z - x k ) #5 (z - x k ) # = # (z - x k ) S (z - x k ) # = 0, e > 0, 

#5 (z - x k ) #{z- x k ) = #5(z- x k ) (z - x k ) # = 0, e < 0. 

The minimal losses, corresponding to the optimal quantum strategy are defined by 
the following formula 

a° = n \z\ 2 +h (q h + x; (n k a + a^) (<ppn k - #) , (6.11) 

where A^, Q k are defined by (6.2), (6.3), and n k , S fe by (5.11), (5.12) with ji = max (0, e) ■ 
Let us also obtain the solution to the corresponding time-continuous optimal control 
problem for the quantum open system, described by the linear stochastic differential 
equations (5.17), (5.18) and the quadratic integral criterion 

n{-.\x(r) | 2 :) +£ (u :| x (t) | 2 : — 2Re$w (t) x (t) + i? 1 \u(t) | 2 ) dt. 

This criterion is obtained by setting oo ~ luA, d ~ i?A, d l ~ in the conditions of 
the Theorem 6.1, and passing to the limit as A — > 0. So, the solution to the quantum 
optimal control problem for the time continuous quantum open system (5.17), (5.18) with 
quantum white noises v, w is defined as the limit of the solution to the discrete problem 
at A — ► 0. 



(6.9) 



The optimal strategy, obtained in this limit, is obviously linear with respect to the 
optimal estimate z(t) of x(t) as in the classical case [20]: u(t) = —X(t)z(t), where 
A (*) = (pSl (t) - tf) /t? 1 , n (r) = n, and ft (t) satisfies the equation: 

-dil (t) /dt + (a + a)n(t)=w-\\ (t) | 2 tf 1 . 

The optimal estimate z (t) is obtained by coherent measurements, corresponding to the 
case p = max (0, e) in the time-continuous Kalman filter, and the minimal mean square 
losses are defined by the integral 

a ° = n \z | 2 +h (n E + ^ (p (t)a + \ (tj) (<i>pn (t) - #) e (t) dt) . 

In particular, when P = r y = e = a + a>0, v — v — a, ft 1 = $ + 9, <d = u, 
we obtain the solution to the optimal control problem for the quantum open oscillator 
matched with the transmission line (2.3) of the wave resistance 7/2 which was considered 
as the motivating example in §2. In this case the equations (5.17), (5.18) are reduced 
to (2.2), (2.3), where the generalized derivatives v (t) = v(dt)/dt, y (t) = y(dt)/dt 
represent the direct and reverse waves on the input of the open oscillator. 

A APPENDIX 

Let A, B be von-Neumann algebras, i.e. selfadjoint weakly closed subalgebras of operators 
in a complex Hilbert space 7i including the identity operator 1, and V, 1Z be predual spaces 
of ultra weakly continuous functionals on A and B, respectively. The elements tt e V 
and p G 1Z are called states on A and B respectively if < n, a > > 0, < p, b > > Va > 
0, b > (a, b > means the non-negative definiteness of the operators a e A and b e B), 
and if < 7r, 1 > = 1, < p, 1 > = 1. Linear operators transforming operators b G B into 
operators a & A are called superoperators, and the predual linear maps V — > TZ are called 
operations. The typical example of a superoperator gives a representation b 1— > u*bu, 
where u is a unitary operator. An operation M : n 1— > 7rM G 72. is called the (statistical) 
morphism if the dual superoperator b 1— > M6 G ^4 is positive 1 M6 > 0, V6 > and 
Ml = 1 (it is convenient to denote the morphisms and dual superoperators by the same 
symbol with the right and the left action respectively: < 7rM, b >=< n, Mb >.) 

A B- valued measure b (d() on some Borel space Z 3 ( is called Z-measurement, if 
b (d() > for any Borel d( C Z and / b (d() = 1 in the same sense. If M : V — > 72 
is a morphism describing a quantum channel, 7r-the state on its input and b [dQ -the 
measurement on its output, then the probability distribution on Z is calculated by any 
of the formulas 

P (d() =< ttM, b {d() >=< tt, Mb (d() > . (A.l) 

Let, for instance, the subalgebras A and B be generated by the operators x and y 
respectively with the canonical commutation relations 

[x,y] = 0, [x,x*]=hl, [y,x*]=>yhl, [y,y*]=ehl, 

1 For a physical realization of the statistical morphisms by conditional expectations of the represen- 
tations a stronger condition of complete positivity [Mbik] i k=1 n > 0,Vn, where [6ifc] i k=1 n > is any 
non-negative definite operator-matrix with the elements bik £ B, should be imposed on the morphisms. 



where 7 G C, e G R and ft > is a constant. 

It may be assumed that y = 7X + v holds, where v is an operator in Ti commuting 
with x and x*, but not commuting with the adjoint one: [v,v*] = (e— | 7 | 2 ) HI, and the 
algebra generated by the pair x, y can be represented in the form of the tensor product 
A <E> B°, where B° is the von-Neumann algebra generated by the operator v. 

We shall write the operators, generated by the operators x and v in the form #ip (x, v) 
where ip (£, 77) are complex-valued functions of £, 77 G C, called symbols, and the notation 
# ■ # indicates such order of action for the operators between them, that first act the 
operators x,v, and then their conjugate. For instance, # | x | 2 # = x*x. In a sufficiently 
wide class of symbols any operator from A® B° can be represented in such a form, and 
this representation is single-valued and injective. In the case y = 7a; + v the operators 
a G A are described by the symbols </?(£, 77) — and the operators 6 G £> by the 

symbols (p (£,77) = (3 (7^ — r/)as in the classical commutative case h = 0. The states in 
this quasi-classical representation are described by distributions (7 (£,77), generalizing the 
probability densities and representing the density operators as the symbols of the con- 
trary order, which are dual to the order for the symbols <p (£,7/). Due to h > 0, x*/Vh 
is the standard creation operator, and x/Vh is the standard annihilation operator, so 
that the representation a = j^a (x) # of operators a G A is normal [19], described by the 
holomorphic symbols a (£) with respect to both £, £. The corresponding symbols of 
the states n on A are described by the Glauber distributions which are defined as 

the linear functionals 

<n,a>=Jp(Qa(Qdi (df = dRe^M) , 

describing the symbols of the density operator 7r, appropriate to the antinormal order. 
The normal order is denoted by the parentheses : : , so we have #a (x) # =: a (x) : 
when [x, x*] > 0. Note, that the antinormal symbol p (£) of the density operator 7r and 
the normal symbol p° (£) are connected by the convolution [21] 

p° (0 = I exp {- I £-£T /n}p (C 1 ) rf^ 1 . (A.2) 

The appropriate representation of the algebra B° ', and hence ^4®i3°, is normal only if 
e >| 7 | 2 , when [t>, t>*] > 0. If m (77) is a distribution which defines a state on £>° and there 
is no statistical dependence, a state on A® B° is described by the product p (£) m (77) and 
a state on the sub-algebra B by the convolution 

r(77) = /m(77- 7 0p(0^- (A.3) 

A superoperator B ^ A, which is dual to a morphism (A.3), is described by the symbol 
transformation 

a (0 = J P (v) m(rj- 7$) ^77. 

For the normality of the appropriate representation b = #/3 (y) # of the operators 
b G B with the distribution (A.3) being Glauber, it is sufficient, that e > 0. When 
e < 0, the distribution r (77) is the normal symbol of the appropriate density operator 
P = #r(y)#. 



Let us consider the complex measurements, described by the measurements of the sum 
Ky + w = z, where w is an operator in TC, which commutes with y and y*, but does not 
commute with the adjoint one w*: 



so that [z, z*\ = (it is assumed that the space H is chosen sufficiently wide, otherwise 
such an operator in 7i may not exist.) 

If i] (() is a distribution describing a state on the algebra B 1 generated by the operator 
w, then the probability distribution of the results of such a measurement on the output of 
the channel is described by the normed with respect to the Lebesgue measure d( density 



In accordance with formula (A.l) such a measurement is described by the B -valued 
measure 



[w, w*] = — e | k | 2 hi, 






(A.4) 




(A.5) 



in accordance with the inequality w*w > {[w*,w] ,0}. When e > and representation 
(A.4) is normal, inequality (A.5) prohibits, in particular, distributions of Dirac 5 -form. 
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